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Abstract 

The spin-orbit interaction generally leads to spin splitting (SS) of electron and hole energy 
states in solids, a splitting that is characterized by a scaling with the wavevector k. Whereas for 
3D bulk zincblende solids the electron (heavy hole) SS exhibits a cubic (linear) scaling with k, in 
2D quantum-wells the electron (heavy hole) SS is currently believed to have a mostly linear (cubic) 
scaling. Such expectations are based on using a small 3D envelope function basis set to describe 
2D physics. By treating instead the 2D system explicitly in a multi-band many-body approach we 
discover a large linear scaling of hole states in 2D. This scaling emerges from hole bands coupling 
that would be unsuspected by the standard model that judges coupling by energy proximity. This 
discovery of a linear Dresselhaus k-scaling for holes in 2D implies a different understanding of 
hole-physics in low-dimensions. 

PACS numbers: 71.70.Ej, 73.21.Fg, 71.15.-m 
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Spin-orbit interaction causes the energy levels of 3D bulk-solids [l| and 2D quantum 
wells (QW's) [2] to exhibit a zero-field SS for sufficiently low-symmetry states. On the 
experimental side, attention has recently focused on spin of holes in 2D quantum-wells 
because of their spin Hall-effect jsj] and in OD quantum dots because of the highly coherent 
hole spin |4j and unusually long hole spin lifetime , promising potential interesting 
applications in spintronic devices and solid state quantum computers. On the theoretical 
side, the long-standing tradition has been to describe hole or electron spin physics 

in low-dimensional (< 3D) nanostructures by an expansion in a rather small basis of 3D 
bulk envelope functions, using effective-mass approaches. In general, when the basis set is 
restricted, the resolution of the expansion is limited. Such low-resolution expansions can be 

n n 

"far sighted" [8j in that the actual atomistic symmetry of the nano object |9| is replaced 
by a fictitious higher symmetry, thus missing important degeneracy-splitting and inter-band 
coupling effects. The farsightedness can be cured by systematically increasing the basis set 
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10] or by introducing ad-hoc terms in the Hamiltonian intended to lower the symmetry 



12]. Both modifications come at the expense of introducing more parameters that are 



not calculable by the envelope function theory itself. Indeed, in the standard model for SS 
of nanostructures 13), 14,[l5|, one uses a phenomenological Hamiltonian where one needs to 
decide at the outset, which 3D bands couple in 2D by the spin-orbit interaction (SOI), rather 
than have the theory force such realization upon us. The potential of missing important 
physical interactions not selected to be present in the model Hamiltonian can be substantial 

The current state of the art for the hole states in 2D is illustrated by the work of Bulaev 
and Loss 13] . Starting from a bulk 3D Hamiltonian restricted to heavy-hole (HH) and light- 



hole (LH) bands ("4 x 4"), they have derived an effective 2x2 Hamiltonian for the 2D heavy 
hole (hhO) subband, demonstrating an exact cancellation of the linear-in-k (Dresselhaus) 
terms [3]. This result (implying a pure, uncoupled hhO state in low-dimensions) has been 
used in numerous theories of hole spin in 2D quantum-well and OD quantum-dot systems, 
including in estimation of hole spin relaxation time 5], demonstration of intrinsic hole spin 
Hall Effect [la ], and other hole spin related phenomena 17, 18, Q, [2^]. We adopt instead 



a different approach in which the 2D nanostructures is viewed as a system in its own right, 
rather than express it in terms of a pre-selected basis drawn from a 3D system. We do so by 
solving the 2D band structure using explicitly the microscopic potential of the 2D system 
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under consideration, thus freeing us from the need to judge at the outset which selected 3D 



bands (e.g., 4 x 4 in Ref. Il3l ) will couple in 2D. The results show that the linear term for holes 
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2JJ. This 



is of the same order of magnitude as the well known linear term for electrons 
discovery of a linear Dresselhaus k-scaling for holes in 2D implies a different understanding 
of hole-physics in low-dimensions. 

The central point of the approach utilized here is that the 3D and 2D systems are each 
described by their own microscopic Hamiltonian which is solved in basis sets whose sole 
property is that it produces a converged solution to the system at hand. Thus, the solution 
reflects only the underlying microscopic Hamiltonian, unmasked by issues of choices of bases 
or pre-selected Hamiltonian terms. We use a rather general microscopic Hamiltonian in the 
U GW representation" 

H GW = -^V 2 + H so + [V ext + V H + V\, (1) 

where V ex t is the electron-ion potential, Vh is the interelectronic Hartree potential of the 
specific (3D or 2D) system; E = iG°W is the self-energy with W being the screened Coulomb 
interaction and G° — 1/ (uj — H°±ie) is the Green's function of a noninteracting Hamiltonian 
H°. The physics represented in Eq. (JTJ includes the atomistic symmetry of the specific (2D 
or 3D) system specified by atomic position vectors in V ext (R, n — d QJ -) and incorporates 
self-consistent electrostatic and exchange-and-correlation effects. This quasi-particle self- 
consistent GW (QSGW) scheme has been established as capable of predicting accurate bulk 



energy bands 



24]. 



22l | . including the Dresselhaus splitting in bulk GaAs |23l . 
The approach described above is computationally intensive and can be readily applied 
only to rather small nanostructures. Thus, for computational expediency, when considering 
larger period quantum wells [e.g, (GaAs) n / (AlAs)„ with n ^> 2] we will map both the small- 
n behavior and the n = oo (bulk) QSGW solutions of Eq. (1) to a screened pseudopotential 
Hamiltonian that captures the former limits yet can be readily applied to orders of magnitude 
larger systems (10 6 atoms were demonstrated in Ref. |8|): 

H PP = -^V 2 + H so + v a (r - R n - d aJ ). (2) 

n,a,j 

Here, the external (V ext ) plus screened (Vh + S) terms of Eq. (1) are described by a su- 
perposition of atom-centered functions v a (where d Q j is the position of atom j of type a in 
the n-th cell R n ). They can be constrained to fit approximately yet accurately the QSGW 



calculated SS of bulk solids (n = oo) 



251 ] and of low-period (n ~ 2) 2D quantum-wells. In 



addition, they reproduce well not only the bulk band gaps throughout the Brillouin zone, 
but also the electron and hole effective-mass tensors, as well as the valence band and con- 
dnction band offsets between the well and barrier tnaterlals QflQ. The spin fitting 
Aj of band % obtained by the direct calculation of Eq. (1) or Eq. (2) will be fitted to the 
conventional form Aj = aik + ^k 3 for i = electrons (e) or holes (h), allowing both a linear 
in k and a cubic in k terms to be present. 

Results of the many-body multi-band calculation: For 3D bulk GaAs, Fig. 1 compares 
our multi-band approach to the ab initio QSGW for 3D bulk GaAs, showing the linear 
and cubic SS of the three lowest conduction bands (CB1, CB2, and CB3) and three highest 
valence bands (HH, LH, and SO). The coefficients of linear and cubic terms as obtained by 
fitting the calculated SS to Aj(/c) are given for each band in the inset to Fig. 1. We find 
that the SS of all bands has a cubic term 7j 3D ' ) , but only HH, LH, and CB3 with angular 
momentum J = 3/2 have linear term af D \ For these bands the cubic and linear splittings 
have opposite signs. The screened pseudopotential solution of Eq. (2) gives similar results 
to QSGW (see inset to Fig. 1). The exception is that Eq. (2) gives ^ = for all bands 



io|. 



because this approach is coreless and this term results from coupling to the core states 
Fortunately, af ^ does not matter much for 2D so we can safely use this method for larger 
systems. 

For 2D (GaAs)2/(AlAs)2 superlattice, the SS of electrons obtained by our atomistic 
multi-band approach agrees well with kp [2l| in that the linear term a^ D ^ oc ai /d 2 
originates from the folding-in of 3D bulk cubic term af^ due to the confinement to a well 
of width d. This is not the case for holes. The SS of 2D valence subbands (hhO, lhO, and 
hhl) is presented in Fig. 2; the first two lines of Table I give the linear and cubic coefficients 
ci^ho an d Ihhti- Both atomistic multi-band methods [Eqs. (1) and (2)] show: (i) a linear 
scaling of SS in addition to a cubic scaling for all three valence subbands including hhO, 
in contrast to only a cubic scaling of hhO in the model-Hamiltonian derived by Bulaev and 
Loss 13]. (ii) The linear term dominates the SS at small k; or^ , cnj^\ and orj^P are 
comparable. 

We next consider a few possible scenarios that might have led to a strong 2D hole splitting, 
finding all but the last to be unlikely. 
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(i) The standard 2x2 model Hamiltonian for 2D does not explain the results. Bulaev 



and Loss 



13j have shown that starting from a 4 x 4 basis in 3D there is no linear terms for 



hhO in the 2D model Hamiltonian. However, they did not include the 3D relativistic cubic 



terms of r 8 „ bands in their derivation. Rashba and Sherman [15| demonstrated earlier that 
the 3D relativistic cubic terms can give rise to 2D linear term for hhO subbands [Eq. (8) 
in Ref. Q. We tested this idea calculating the 2D SS using a pseudopotential which was 
constructed to fit the Rashba-Sherman 3D band structure including the relativistic cubic 
terms. The results of the pseudopotential calculation for 2D (in which all bands are allowed 
to couple) are compared with the Rashba-Sherman 2D model (in which the 2D hhO band is 
uncoupled). Fig. 3 and Table I show that for sufficiently large superlattice periods n, for 



which the 2D model of Ref. 



la is applicable, the model recovers only a small fraction of the 



multi-band results for linear SS of 2D hhO. For smaller periods (n < 20) we report in Fig. 3 
a clear non-monotonic period dependence of both a^ho and ctf^\ with a^ho to ct~^ ratios 
varying from ~ 10 to < 1 as the period decreases from 50 to 2 ML. This is in sharp contrast 
with the predictions of the model Hamiltonian which predicts a monotonic increase of 
linear terms and a ratio of ot^M ^° a< iho^ which is independent of d. 

(ii) Interfaces induce only minor 2D linear splitting: In 2D quantum- wells, Foreman [U 
suggested an interface induced linear term that originates from valence band coupling to the 
p-like T 4 states. For 2D GaAs/AlAs quantum well, he 12j estimated that independent of 
period this interfacial linear term is in the range of 20 — 30 meVA. In contrast, we find for 
near-bulk (large-period) superlattice such as (GaAs)so/(AlAs)2o (Table I) that aj^Q is very 
close to its bulk value of QSGW. Thus, even if we assume that the interfacial linear term 
is solely responsible for the remai ning discrepancy from bulk value, it is even smaller than 



what was estimated by Foreman 12j. We conclude that the interfacial linear term, if it 
exists, is only a minor contribution to ct^m ■ 

(Hi) Core-valence coupling is not the reason either: One might have suspected that 
quantum confinement pushes the valence bands closer to core levels and thus increases 
the linear term due to increased coupling. This too can be excluded since our screened 
pseudopotential is a coreless method, yet we find similar results as the all-electron QSGW 
calculation. 

(iv) Undiscovered spin-orbit linear terms: The quantitative and qualitative disagree- 
ments of atomistic multi-band calculation with the standard model Hamiltonian approach 
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suggest possible undiscovered SOI terms, which are not included in the model Hamiltonian. 
Such terms due to symmetry lowering down from 3D T<j bulk symmetry to 2D Z^d quantum 
well symmetry could originate from coupling of 3D bands via the 2D potential and SOI. Such 
coupling is signaled by the distinctly nonparabolic 2D energy dispersion curves manifesting 



clear anti-crossings between neighboring subbands, mu 
inset to Fig. 3. In the model Hamiltonian approach 13 



ti-band calculation displayed in the 



15] the hhO, hhl, wavefunctions 



near zone-center in 2D all derive from a single bulk state |HH) and similarly, all lhO, lhl 
wavefunctions in 2D derive from a single bulk state |LH). However, the inset to Fig. 3 shows 
that the lhO level lies between the hhO and hhl levels in 2D for all multi-band calculated 
GaAs periods. Clearly, the coupling of 3D states and its effects on SS of 2D bands can not 



be ignored as done in the model Hamiltonian approach for extremely small-k range 
A better approach than this "band decoupled" model Hamiltonian 



13 



13 



15| 



15J approximation 



allows the 2D state hhO to derive from a few bulk states. In such a "mixing of decoupled 
states" approximation, 

= «ffi(HH)|HH) + «©(LH)|LH) + • • • , (3) 

where iujj^(A) is the percent weig ht of 3D state A = HH, LH, ... in the 2D state 
In the band decoupled model wj^\\) = for A ^ HH. We have calculated the weights 
by numerical projection of the 2D pseudopotential wavefunctions and show them in Fig. 4. 
In striking contrast to the assumption made in model Hamiltonian approach, the interband 
coupling is large even at zone center. We see that for long-period (GaAs thickness ^ 20 ML) 
^hhn * s ma de of 90% |HH) and 5% |LH), but for shorter periods the mixing increases: the 
HH content drops to ~ 70 — 80 % and the LH content raises to 10 — 20 %. This monotonicly 
enhanced mixing of LH and HH into the 2D hhO as the QW period is reduced signals the 
breakdown of the model Hamiltonian thinking that neglects such mixing on the ground 
that the energy splitting of hhO-lhO must be larger than that of hhO-hhl for sufficient small 
periods [kJ. 

The linear coefficient of 2D hhO SS can be written in terms of the weights in Eq. (3) in 
a model of "mixing of decoupled states" as 

J?D) _ (2D) mm ~ (2D) (2D) (J m ~(2D) m 

a hho - w hho \ nn ) a hho + w hho \ Lti ) a iho H > l 4 J 

where &hhn(®iho ) * s ^ ne contribution of a single bulk HH (or LH) band to linear SS of 



2D hhO subband, which had been derived by Rashba and Sherman 



151 ] (the negative sign 



accounts for band repulsion effect). The result of the first two terms in Eq. (4) is shown as 
open squares in Fig. 3 and are compared with the Muti-band calculated (solid circles). 
We see that the mixing of decoupled states [Eq. (4)] gives a much better approximation to 
the full calculation than the model Hamiltonian treating one decoupled band at the time 
(open circles in Fig.3) 26(|. Thus, the mixing of bulk bands leads to a large linear SS of 2D 
hhO, and is unsuspected by the standard model that judge coupling by energy proximity. 

The emergence of a large linear term for Dresselhaus hole SS in 2D nanostructures sug- 
gests (i) the dominance of Dresselhaus over Rashba SOI (having a cubic term as its lowest 
order term) 13|. (ii) a larger spin-Hall effect jla ]. and (iii) an explanation of the observed 
,arge optica! anisotropy Q. The occurrence of a larger SS of hhO correspond™ to HH- 
LH coupling leads to a short hole spin-relaxtion time in 2D quantum-wells |6|] from the 
D'Yakonov and Perel (DP) mechanism [zj]. 
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was funded by the U.S. Department of Energy, Office of Science, Basic Energy Science, 
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MvS was supported by ONR, project N00014-07-1-0479. 
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TABLE I: The coefficients of linear and cubic terms and J^ 2D ^ for 2D superlattice as well as 

3D bulk GaAs. Here a\ 2D ^ units in meVA and "f[ 2D ^ in eVA 3 . The symbol in parentheses denotes 
the character of the folded-in band (bulk T-like or X-like) at the Brillouin-zone center F of the 
2D system. We show results from QSGW [Eq. (1)] as well as from screened pseudopotential [Eq. 
(2)]. The pseudopotential used for 2D system in this Table was constructed specifically to predict 
comparable bulk SS to model Hamiltonian result l^] (see text) rather than QSGW values. 



Superlattice 
(GaAs) n /(AlAs) m 


CBl 

(2D) (2D) 


VBl 

(2D) (2D) 
a hh0 ^hhO 


GW: 2/2 (X) 
PP- 2/2 (X) 

.1.1. > _■ j _ ■ I ... \_ } 


61.5 0.2 


102.6 151.3 
33Q 6 15 4 


PP: 2/4 (X) 
PP: 2/10 (X) 
PP: 2/20 (X) 


8.5 0.2 
23.0 2.2 
20.0 10.5 


352.9 1.5 

379.1 13.2 

420.2 13.2 


PP: 20/20 (r) 
PP: 30/20 (r) 
PP: 50/20 (r) 
PP: 80/20 (r) 


71.7 12.2 
41.5 7.3 
17.9 3.6 
7.7 1.5 


206.1 95.9 
114.5 111.8 
50.4 126.6 
22.1 113.6 


PP: 3D GaAs 
GW: 3D GaAs 
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0.0 8.3 
12.6 3.1 
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FIG. 1: (color online) QSGW predicted SS of 3D bulk GaAs for the three lowest conduction 
bands (a) CB1, (b) CB2, and (c) CB3, and three highest valence bands (d) VB1, (e) VB2, and 
(f) VB3. Dresselhaus constants ctj in meVA and ji in eVA 3 predicted by QSGW and screened 
pseudopotential fit to QSGW, respectively, are given for each band in the inset. 
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FIG. 2: QSGW [Eq. (1)] predicted SS (red line) of 2D (GaAs) 2 /(AlAs) 2 superlattice for the 
three highest valence subbands (a) hhO, (b) lhO, and (c) hhl. The corresponding SS calculated by 
screened pseudopotential fit to model Hamiltonian result [Eq. (2)] is also given in black line 
for each subband. 
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FIG. 3: (a) Comparison of linear SS coefficients a hh0 and a lh0 calculated from present pseu- 
dopotential multi-band approach (same as in Table I; solid circles), as well as by the Rashba- 
Sherman model Hamiltonian (open circles), and mixing approximation [Eq. (4); open squares] for 
2D (GaAs) n /(AlAs)2o quantum wells, (b) Expanded scale for long-period QW's highlighting the 
comparison of hhO subbands for the direct calculation vs model Hamiltonian and mixing approxi- 
mation. Inset shows energy dispersion of valence subbands for 2D (GaAs) n /(AlAs)2o calculated by 
atomistic pseudopotential approach. The two spin subbands for each orbit subband are represented 
by low energy red line and a high energy black line. 
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FIG. 4: Percent of bulk HH and LH bands mixed into the 2D hhO and lhO subbands at zone center 
calculated by numerical projection. The results are almost same in a small in-plane k range [e.g., 
k < 0.01(27r/a)]. Note that the sum of bulk HH and LH content is less than 100%. In our atomistic 
multi-band calculation, the subband ordering from top to bottom is hhO, lhO, hhl, ... for all GaAs 
period. 
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